Seasonal Patterns of Mpox Index Cases, Africa, 1970-2021.

Across 133 confirmed mpox zoonotic index cases reported during 1970-2021 in Africa, cases occurred year-round near the equator, where climate is consistent. However, in tropical regions of the northern hemisphere under a dry/wet season cycle, cases occurred seasonally. Our findings further support the seasonality of mpox zoonotic transmission risk.

Identifying the Most Probable Mammal Reservoir Hosts for Monkeypox Virus Based on Ecological Niche Comparisons.

Previous human cases or epidemics have suggested that Monkeypox virus (MPXV) can be transmitted through contact with animals of African rainforests. Although MPXV has been identified in many mammal species, most are likely secondary hosts, and the reservoir host has yet to be discovered. In this study, we provide the full list of African mammal genera (and species) in which MPXV was previously detected, and predict the geographic distributions of all species of these genera based on museum specimens and an ecological niche modelling (ENM) method. Then, we reconstruct the ecological niche of MPXV using georeferenced data on animal MPXV sequences and human index cases, and conduct overlap analyses with the ecological niches inferred for 99 mammal species, in order to identify the most probable animal reservoir. Our results show that the MPXV niche covers three African rainforests: the Congo Basin, and Upper and Lower Guinean forests. The four mammal species showing the best niche overlap with MPXV are all arboreal rodents, including three squirrels: Funisciurus anerythrus, Funisciurus pyrropus, Heliosciurus rufobrachium, and Graphiurus lorraineus. We conclude that the most probable MPXV reservoir is F. anerythrus based on two niche overlap metrics, the areas of higher probabilities of occurrence, and available data on MPXV detection.

Detection and genetic characterization of circoviruses in more than 80 bat species from eight countries on four continents.

Several bat-associated circoviruses and circular rep-encoding single-stranded DNA (CRESS DNA) viruses have been described, but the exact diversity and host species of these viruses are often unknown. Our goal was to describe the diversity of bat-associated circoviruses and cirliviruses, thus, 424 bat samples from more than 80 species were collected on four continents. The samples were screened for circoviruses using PCR and the resulting amino acid sequences were subjected to phylogenetic analysis. The majority of bat strains were classified in the genus Circovirus and some strains in the genus Cyclovirus and the clades CRESS1 and CRESS3. Some strains, however, could only be classified at the taxonomic level of the order and were not classified in any of the accepted or proposed clades. In the family Circoviridae, 71 new species have been predicted. This screening of bat samples revealed a great diversity of circoviruses and cirliviruses. These studies underline the importance of the discovery and description of new cirliviruses and the need to establish new species and families in the order Cirlivirales.

Retracing Phylogenetic, Host and Geographic Origins of Coronaviruses with Coloured Genomic Bootstrap Barcodes: SARS-CoV and SARS-CoV-2 as Case Studies.

Phylogenetic trees of coronaviruses are difficult to interpret because they undergo frequent genomic recombination. Here, we propose a new method, coloured genomic bootstrap (CGB) barcodes, to highlight the polyphyletic origins of human sarbecoviruses and understand their host and geographic origins. The results indicate that SARS-CoV and SARS-CoV-2 contain genomic regions of mixed ancestry originating from horseshoe bat (Rhinolophus) viruses. First, different regions of SARS-CoV share exclusive ancestry with five Rhinolophus viruses from Southwest China (RfYNLF/31C: 17.9%; RpF46: 3.3%; RspSC2018: 2.0%; Rpe3: 1.3%; RaLYRa11: 1.0%) and 97% of its genome can be related to bat viruses from Yunnan (China), supporting its emergence in the Rhinolophus species of this province. Second, different regions of SARS-CoV-2 share exclusive ancestry with eight Rhinolophus viruses from Yunnan (RpYN06: 5.8%; RaTG13: 4.8%; RmYN02: 3.8%), Laos (RpBANAL103: 3.3%; RmarBANAL236: 1.7%; RmBANAL52: 1.0%; RmBANAL247: 0.7%), and Cambodia (RshSTT200: 2.3%), and 98% of its genome can be related to bat viruses from northern Laos and Yunnan, supporting its emergence in the Rhinolophus species of this region. Although CGB barcodes are very useful in retracing the origins of human sarbecoviruses, further investigations are needed to better take into account the diversity of coronaviruses in bats from Cambodia, Laos, Myanmar, Thailand and Vietnam.

Phylogeny of sea spiders (Arthropoda: Pycnogonida) inferred from mitochondrial genome and 18S ribosomal RNA gene sequences.

The phylogeny of sea spiders has been debated for more than a century. Despite several molecular studies in the last twenty years, interfamilial relationships remain uncertain. In the present study, relationships within Pycnogonida are examined in the light of a new dataset composed of 160 mitochondrial genomes (including 152 new sequences) and 130 18S rRNA gene sequences (including 120 new sequences), from 141 sea spider morphospecies representing 26 genera and 9 families. Node congruence between mitochondrial and nuclear markers was analysed to identify the most reliable relationships. We also reanalysed a multilocus dataset previously published and showed that the high percentages of missing data make phylogenetic conclusions difficult and uncertain. Our results support the monophyly of most families currently accepted, except Callipallenidae and Nymphonidae, the monophyly of the superfamilies Ammotheoidea (Ammotheidae + Pallenopsidae), Nymphonoidea (Nymphonidae + Callipallenidae), Phoxichilidioidea (Phoxichilidiidae + Endeidae) and Colossendeoidea (Colossendeidae + Pycnogonidae + Rhynchothoracidae), and the sister-group relationship between Ammotheoidea and Phoxichilidioidea. We discuss the morphological evolution of sea spiders, identifying homoplastic characters and possible synapomorphies. We also discuss the palaeontological and phylogenetic arguments supporting either a radiation of sea spiders prior to Jurassic or a progressive diversification from Ordovician or Cambrian.

National Monkeypox Surveillance, Central African Republic, 2001-2021.

We analyzed monkeypox disease surveillance in Central African Republic (CAR) during 2001-2021. Surveillance data show 95 suspected outbreaks, 40 of which were confirmed as monkeypox, comprising 99 confirmed and 61 suspected monkeypox cases. After 2018, CAR's annual rate of confirmed outbreaks increased, and 65% of outbreaks occurred in 2 forested regions bordering the Democratic Republic of the Congo. The median patient age for confirmed cases was 15.5 years. The overall case-fatality ratio was 7.5% (12/160) for confirmed and suspected cases, 9.6% (8/83) for children <16 years of age. Decreasing cross-protective immunity from smallpox vaccination and recent ecologic alterations likely contribute to increased monkeypox outbreaks in Central Africa. High fatality rates associated with monkeypox virus clade I also are a local and international concern. Ongoing investigations of zoonotic sources and environmental changes that increase human exposure could inform practices to prevent monkeypox expansion into local communities and beyond endemic areas.

Draft genome of the lowland anoa (Bubalus depressicornis) and comparison with buffalo genome assemblies (Bovidae, Bubalina).

Genomic data for wild species of the genus Bubalus (Asian buffaloes) are still lacking while several whole genomes are currently available for domestic water buffaloes. To address this, we sequenced the genome of a wild endangered dwarf buffalo, the lowland anoa (Bubalus depressicornis), produced a draft genome assembly and made comparison to published buffalo genomes. The lowland anoa genome assembly was 2.56 Gbp long and contained 103,135 contigs, the longest contig being 337.39 kbp long. N50 and L50 values were 38.73 and 19.83 kbp, respectively, mean coverage was 44× and GC content was 41.74%. Two strategies were adopted to evaluate genome completeness: (1) determination of genomic features with de novo and homology-based predictions using annotations of chromosome-level genome assembly of the river buffalo and (2) employment of benchmarking against universal single-copy orthologs (BUSCO). Homology-based predictions identified 94.51% complete and 3.65% partial genomic features. De novo gene predictions identified 32,393 genes, representing 97.14% of the reference's annotated genes, whilst BUSCO search against the mammalian orthologs database identified 71.1% complete, 11.7% fragmented, and 17.2% missing orthologs, indicating a good level of completeness for downstream analyses. Repeat analyses indicated that the lowland anoa genome contains 42.12% of repetitive regions. The genome assembly of the lowland anoa is expected to contribute to comparative genome analyses among bovid species.

Variation in synonymous nucleotide composition among genomes of sarbecoviruses and consequences for the origin of COVID-19.

The subgenus Sarbecovirus includes two human viruses, SARS-CoV and SARS-CoV-2, respectively responsible for the SARS epidemic and COVID-19 pandemic, as well as many bat viruses and two pangolin viruses. Here, the synonymous nucleotide composition (SNC) of Sarbecovirus genomes was analysed by examining third codon-positions, dinucleotides, and degenerate codons. The results show evidence for the eight following groups: (i) SARS-CoV related coronaviruses (SCoVrC including many bat viruses from China), (ii) SARS-CoV-2 related coronaviruses (SCoV2rC; including five bat viruses from Cambodia, Thailand and Yunnan), (iii) pangolin sarbecoviruses, (iv) three bat sarbecoviruses showing evidence of recombination between SCoVrC and SCoV2rC genomes, (v) two highly divergent bat sarbecoviruses from Yunnan, (vi) the bat sarbecovirus from Japan, (vii) the bat sarbecovirus from Bulgaria, and (viii) the bat sarbecovirus from Kenya. All these groups can be diagnosed by specific nucleotide compositional features except the one concerned by recombination between SCoVrC and SCoV2rC. In particular, SCoV2rC genomes have less cytosines and more uracils at third codon-positions than other sarbecoviruses, whereas the genomes of pangolin sarbecoviruses show more adenines at third codon-positions. I suggest that taxonomic differences in the imbalanced nucleotide pools available in host cells during viral replication can explain the eight groups of SNC here detected among Sarbecovirus genomes. A related effect due to hibernating bats and their latitudinal distribution is also discussed. I conclude that the two independent host switches from Rhinolophus bats to pangolins resulted in convergent mutational constraints and that SARS-CoV-2 emerged directly from a horseshoe bat sarbecovirus.

Presence of Recombinant Bat Coronavirus GCCDC1 in Cambodian Bats.

Bats have been recognized as an exceptional viral reservoir, especially for coronaviruses. At least three bat zoonotic coronaviruses (SARS-CoV, MERS-CoV and SARS-CoV-2) have been shown to cause severe diseases in humans and it is expected more will emerge. One of the major features of CoVs is that they are all highly prone to recombination. An extreme example is the insertion of the P10 gene from reoviruses in the bat CoV GCCDC1, first discovered in Rousettus leschenaultii bats in China. Here, we report the detection of GCCDC1 in four different bat species (Eonycteris spelaea, Cynopterus sphinx, Rhinolophus shameli and Rousettus sp.) in Cambodia. This finding demonstrates a much broader geographic and bat species range for this virus and indicates common cross-species transmission. Interestingly, one of the bat samples showed a co-infection with an Alpha CoV most closely related to RsYN14, a virus recently discovered in the same genus (Rhinolophus) of bat in Yunnan, China, 2020. Taken together, our latest findings highlight the need to conduct active surveillance in bats to assess the risk of emerging CoVs, especially in Southeast Asia.

Genomic Bootstrap Barcodes and Their Application to Study the Evolution of Sarbecoviruses.

Recombination creates mosaic genomes containing regions with mixed ancestry, and the accumulation of such events over time can complicate greatly many aspects of evolutionary inference. Here, we developed a sliding window bootstrap (SWB) method to generate genomic bootstrap (GB) barcodes to highlight the regions supporting phylogenetic relationships. The method was applied to an alignment of 56 sarbecoviruses, including SARS-CoV and SARS-CoV-2, responsible for the SARS epidemic and COVID-19 pandemic, respectively. The SWB analyses were also used to construct a consensus tree showing the most reliable relationships and better interpret hidden phylogenetic signals. Our results revealed that most relationships were supported by just a few genomic regions and confirmed that three divergent lineages could be found in bats from Yunnan: SCoVrC, which groups SARS-CoV related coronaviruses from China; SCoV2rC, which includes SARS-CoV-2 related coronaviruses from Southeast Asia and Yunnan; and YunSar, which contains a few highly divergent viruses recently described in Yunnan. The GB barcodes showed evidence for ancient recombination between SCoV2rC and YunSar genomes, as well as more recent recombination events between SCoVrC and SCoV2rC genomes. The recombination and phylogeographic patterns suggest a strong host-dependent selection of the viral RNA-dependent RNA polymerase. In addition, SARS-CoV-2 appears as a mosaic genome composed of regions sharing recent ancestry with three bat SCoV2rCs from Yunnan (RmYN02, RpYN06, and RaTG13) or related to more ancient ancestors in bats from Yunnan and Southeast Asia. Finally, our results suggest that viral circular RNAs may be key molecules for the mechanism of recombination.

A novel SARS-CoV-2 related coronavirus in bats from Cambodia.

Knowledge of the origin and reservoir of the coronavirus responsible for the ongoing COVID-19 pandemic is still fragmentary. To date, the closest relatives to SARS-CoV-2 have been detected in Rhinolophus bats sampled in the Yunnan province, China. Here we describe the identification of SARS-CoV-2 related coronaviruses in two Rhinolophus shameli bats sampled in Cambodia in 2010. Metagenomic sequencing identifies nearly identical viruses sharing 92.6% nucleotide identity with SARS-CoV-2. Most genomic regions are closely related to SARS-CoV-2, with the exception of a region of the spike, which is not compatible with human ACE2-mediated entry. The discovery of these viruses in a bat species not found in China indicates that SARS-CoV-2 related viruses have a much wider geographic distribution than previously reported, and suggests that Southeast Asia represents a key area to consider for future surveillance for coronaviruses.

Inferring the ecological niche of bat viruses closely related to SARS-CoV-2 using phylogeographic analyses of Rhinolophus species.

The Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) is the causal agent of the coronavirus disease 2019 (COVID-19) pandemic. To date, viruses closely related to SARS-CoV-2 have been reported in four bat species: Rhinolophus acuminatus, Rhinolophus affinis, Rhinolophus malayanus, and Rhinolophus shameli. Here, we analysed 343 sequences of the mitochondrial cytochrome c oxidase subunit 1 gene (CO1) from georeferenced bats of the four Rhinolophus species identified as reservoirs of viruses closely related to SARS-CoV-2. Haplotype networks were constructed in order to investigate patterns of genetic diversity among bat populations of Southeast Asia and China. No strong geographic structure was found for the four Rhinolophus species, suggesting high dispersal capacity. The ecological niche of bat viruses closely related to SARS-CoV-2 was predicted using the four localities in which bat viruses were recently discovered and the localities where bats showed the same CO1 haplotypes than virus-positive bats. The ecological niche of bat viruses related to SARS-CoV was deduced from the localities where bat viruses were previously detected. The results show that the ecological niche of bat viruses related to SARS-CoV2 includes several regions of mainland Southeast Asia whereas the ecological niche of bat viruses related to SARS-CoV is mainly restricted to China. In agreement with these results, human populations in Laos, Vietnam, Cambodia, and Thailand appear to be much less affected by the COVID-19 pandemic than other countries of Southeast Asia. In the climatic transitional zone between the two ecological niches (southern Yunnan, northern Laos, northern Vietnam), genomic recombination between highly divergent viruses is more likely to occur. Considering the limited data and the risk of recombinant bat-CoVs emergence as the source of new pandemics in humans, the bat populations in these regions should be under surveillance.

Genomic history of human monkey pox infections in the Central African Republic between 2001 and 2018.

Monkeypox is an emerging infectious disease, which has a clinical presentation similar to smallpox. In the two past decades, Central Africa has seen an increase in the frequency of cases, with many monkeypox virus (MPXV) isolates detected in the Democratic Republic of Congo (DRC) and the Central African Republic (CAR). To date, no complete MPXV viral genome has been published from the human cases identified in the CAR. The objective of this study was to sequence the full genome of 10 MPXV isolates collected during the CAR epidemics between 2001 and 2018 in order to determine their phylogenetic relationships among MPXV lineages previously described in Central Africa and West Africa. Our phylogenetic results indicate that the 10 CAR isolates belong to three lineages closely related to those found in DRC. The phylogenetic pattern shows that all of them emerged in the rainforest block of the Congo Basin. Since most human index cases in CAR occurred at the northern edge of western and eastern rainforests, transmissions from wild animals living in the rainforest is the most probable hypothesis. In addition, molecular dating estimates suggest that periods of intense political instability resulting in population movements within the country often associated also with increased poverty may have led to more frequent contact with host wild animals. The CAR socio-economic situation, armed conflicts and ecological disturbances will likely incite populations to interact more and more with wild animals and thus increase the risk of zoonotic spillover.

The genome of the lowland anoa (Bubalus depressicornis) illuminates the origin of river and swamp buffalo.

Two types of domestic water buffalo are currently recognized: the river buffalo from the Indian subcontinent and Mediterranean countries and the swamp buffalo from China and Southeast Asia. To test the hypothesis of two separate species of water buffalo, we sequenced the genome of the lowland anoa, Bubalus depressicornis, which is a dwarf wild buffalo endemic to Sulawesi, and two genomes of swamp buffalo, and made comparisons with 12 additional genomes. Three genomic data sets were constructed to infer phylogenetic relationships: the mitochondrial genome (15,468 bp; maternal transmission), two concatenated Y-chromosomal genes, AMELY and DDX3Y (20,036 bp; paternal transmission), and a selection of 30 nuclear genes representing all cattle chromosomes (364,887 bp; biparental transmission). The comparisons between our 30 nuclear gene sequences obtained by read mapping and those directly extracted from Bos taurus and Bubalus bubalis genome assemblies show that the mapping approach revealed higher levels of heterozygosity at both nucleotide sites and indels (insertions and deletions) (0.09-0.15%), as well as several sequence errors (0.07%). Our phylogenetic and molecular dating analyses provide strong evidence that the lowland anoa, river buffalo, and swamp buffalo are three distinct taxa which separated rapidly from each other during the Pleistocene epoch. We therefore conclude that two species of domestic water buffalo should be distinguished: Bubalus bubalis for the river buffalo and Bubalus kerabau for the swamp buffalo. The new classification can have deep implications for understanding the evolution and selection of domesticated forms and for the conservation and management of wild buffalo populations in South and Southeast Asia.

Correction: Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes.

[This corrects the article DOI: 10.1371/journal.pone.0240770.].

Haplotype diversity in the mitochondrial genome of the Egyptian river buffalo (Bubalus bubalis).

Here, we sequenced the complete mitochondrial genome of 29 Egyptian river buffaloes collected in two breeding stations of Egypt. The genome is 16,357-16,359 base pairs in length and contains the 37 genes found in a typical mammalian genome. The overall base composition is A: 33.1%, C: 26.6%, G: 13.9%, and T: 26.4%. Our analyses confirm that the mitochondrial genomes of swamp and river buffaloes are divergent (mean nucleotide distance = 2.3%), and show that Indian river buffalo haplotypes cluster into three haplogroups, named RB1, RB2, and RB3 (mean distance = 0.25-0.26%) and that the 24 Egyptian buffalo haplotypes fall into RB1 (with the Bangladeshi, Chinese and Italian buffalo haplotypes) and RB2.

Evolutionary history of Carnivora (Mammalia, Laurasiatheria) inferred from mitochondrial genomes.

The order Carnivora, which currently includes 296 species classified into 16 families, is distributed across all continents. The phylogeny and the timing of diversification of members of the order are still a matter of debate. Here, complete mitochondrial genomes were analysed to reconstruct the phylogenetic relationships and to estimate divergence times among species of Carnivora. We assembled 51 new mitogenomes from 13 families, and aligned them with available mitogenomes by selecting only those showing more than 1% of nucleotide divergence and excluding those suspected to be of low-quality or from misidentified taxa. Our final alignment included 220 taxa representing 2,442 mitogenomes. Our analyses led to a robust resolution of suprafamilial and intrafamilial relationships. We identified 21 fossil calibration points to estimate a molecular timescale for carnivorans. According to our divergence time estimates, crown carnivorans appeared during or just after the Early Eocene Climatic Optimum; all major groups of Caniformia (Cynoidea/Arctoidea; Ursidae; Musteloidea/Pinnipedia) diverged from each other during the Eocene, while all major groups of Feliformia (Nandiniidae; Feloidea; Viverroidea) diversified more recently during the Oligocene, with a basal divergence of Nandinia at the Eocene/Oligocene transition; intrafamilial divergences occurred during the Miocene, except for the Procyonidae, as Potos separated from other genera during the Oligocene.

A comparative approach for species delimitation based on multiple methods of multi-locus DNA sequence analysis: A case study of the genus Giraffa (Mammalia, Cetartiodactyla).

Molecular data are now commonly used in taxonomy for delimiting cryptic species. In the case of giraffes, which were treated as a single species (Giraffa camelopardalis) during half of a century, several molecular studies have suggested a splitting into four to seven species, but the criteria applied for taxonomic delimitation were not fully described. In this study, we have analysed all multi-locus DNA sequences available for giraffes using multispecies coalescent (MSC: *BEAST, BPP and STACEY), population genetic (STRUCTURE, allelic networks, haplotype network and bootstrapping, haplowebs and conspecificity matrix) and phylogenetic (MrBayes, PhyML, SuperTRI) methods to identify the number of species. Our results show that depending on the method chosen, different taxonomic hypotheses, recognizing from two to six species, can be considered for the genus Giraffa. Our results confirm that MSC methods can lead to taxonomic over-splitting, as they delimit geographic structure rather than species. The 3-species hypothesis, which recognizes G. camelopardalis sensu strico A, G. giraffa, and G. tippelskirchi, is highly supported by phylogenetic analyses and also corroborated by most population genetic and MSC analyses. The three species show high levels of nucleotide divergence in both nuclear (0.35-0.51%) and mitochondrial sequences (3-4%), and they are characterised by 7 to 12 exclusive synapomorphies (ES) detected in nine of the 21 nuclear introns analysed for this study. By contrast, other putative species, such as G. peralta, G. reticulata, G. thornicrofti or G. tippelskirchi sensu stricto, do not exhibit any ES in the nuclear genes. A robust mito-nuclear conflict was found for the position and monophyly of G. giraffa and G. tippelskirchi, which is interpreted as the result of a mitochondrial introgression from Masai to southeastern giraffe during the Pleistocene and nuclear gene flow mediated by male dispersal between southern populations (subspecies G. g. giraffa and G. g. angolensis).

Multi-locus phylogeny of the tribe Tragelaphini (Mammalia, Bovidae) and species delimitation in bushbuck: Evidence for chromosomal speciation mediated by interspecific hybridization.

The bushbuck is the most widespread bovid species in Africa. Previous mitochondrial studies have revealed a polyphyletic pattern suggesting the possible existence of two distinct species. To assess this issue, we have sequenced 16 nuclear genes and one mitochondrial fragment (cytochrome b gene + control region) for most species of the tribe Tragelaphini, including seven bushbuck individuals belonging to the two divergent mtDNA haplogroups, Scriptus and Sylvaticus. Our phylogenetic analyses show that the Scriptus lineage is a sister-group of Sylvaticus in the nuclear tree, whereas it is related to Tragelaphus angasii in the mitochondrial tree. This mito-nuclear discordance indicates that the mitochondrial genome of Scriptus was acquired by introgression after one or several past events of hybridization between bushbuck and an extinct species closely related to T. angasii. The division into two bushbuck species is supported by the analyses of nuclear markers and by the karyotype here described for T. scriptus (2n = 57 M/58F), which is strikingly distinct from the one previously found for T. sylvaticus (2n = 33 M/34F). Molecular dating estimates suggest that the two species separated during the Early Pleistocene after an event of interspecific hybridization, which may have mediated massive chromosomal rearrangements in the common ancestor of T. scriptus.